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ABSTRACT 

Context. It is hypothesized that low-mass young stellar objects undergo eruptive phases during their early evolution. These eruptions 
are thought to be caused by highly increased mass accretion from the disk onto the star, and therefore play an important role in the 
early evolution of Sun-like stars, of their circumstellar disks (structure, dust composition), and in the formation of their planetary 
systems. The outburst of VI647 Ori between 2003 and 2006 offered a rare opportunity to investigate such an accretion event. 
Aims. By means of our interferometry observing campaign during this outburst, supplemented by other observations, we investigate 
the temporal evolution of the inner circumstellar structure of VI647 Ori, the region where Earth-like planets could be bom. We also 
study the role of the changing extinction in the brightening of the object and separate it from the accretional brightening. 
Methods. We observed VI 647 Ori with MIDI on the VLTI at two epochs in this outburst. First, during the slowly fading plateau 
phase (2005 March) and second, just before the rapid fading of the object (2005 September), which ended the outburst. We used the 
radiative transfer code MC3D to fit the interferometry data and the spectral energy distributions from five different epochs at different 
stages of the outburst. The comparison of these models allowed us to trace structural changes in the system on AU-scales. We also 
considered qualitative alternatives for the interpretation of our data. 

Results. We found that the disk and the envelope are similar to those of non-eruptive young stars and that the accretion rate varied 
during the outburst. We also found evidence for the increase of the inner radii of the circumstellar disk and envelope at the beginning 
of the outburst. Furthermore, the change of the interferometric visibilities indicates structural changes in the circumstellar material. 
We test a few scenarios to interpret these data. We also speculate that the changes are caused by the fading of the central source, which 
is not immediately followed by the fading of the outer regions. 

Conclusions. We found that most of our results fit in the canonical picture of young eruptive stars. Our study provided dynamical 
information from the regions of the innermost few AU of the system: changes of the inner radii of the disk and envelope. However, if 
the delay in the fading of the disk is responsible for the changes seen in the MIDI data, the effect should be confirmed by dynamical 
modeling. 

Key words, stars: formation - stars: circumstellar matter - stars: individual: V1647 Ori - infrared: stars - techniques: interferometric 



1. Introduction 

The origin of the matter of low-mass stars is interstellar gas and 
dust, which fall onto the protostar through a circumstellar ac- 
cretion disk. It is becoming increasingly evident that the accu- 
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mulation of stellar material is an episodic process: the typically 
low accretion rate suddenly increases for limited periods, de- 
livering a significant amount of mass onto the stellar surface 
(e.g., Evans et al. 2009) . Many young stars exhibit temporal 
brightenings possibly caused by fluctuating accretion (Herbst et 
al. .1994 , Sicilia-Aguilar et al. 2008 ). The eruptions of the FU 
Orionis- and EX Lupi-type classes of variable stars (FUors and 
EXors) may represent the most intense bursts (e.g., Hartmann 
& Kenvon [T996l hereafter HK96, Herbig 2007D . Such outbursts 
are characterized by an optical brightening of 2-6 magnitude. 
FUor eruptions can last from some years to several decades, 
while EXor-outbursts are usually considered as shorter coun- 
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terparts of those of FUors (i.e., some months long). A num- 
ber of possible scenarios for the triggering mechanism of the 
outburst have been proposed: (a) interactions of binary or mul- 
tiple systems where tidal forces disturb the circumstellar disk 
(Bonnell & Bastien 119921 ). (b) thermal instabilities in the disk 
(Bell et al. 1 19951 1. (c) planet-disk interactions, where thermal in- 
stabilities in the disk are caused by the presence of a massive 
planet (Lodato & Clarke I2004I I. (d) gravitational instabilities in 
the disk (Vorobyov & Basu |25^ Zhu et al. |2009] l. or (e) in- 
teraction between the disk and the magnetic field (D'Angelo & 
Spruit |20T0] |. The difference between the physics of FUor and 
EXor-outbursts is unclear. Although currently only about 20 ob- 
jects are classified as eruptive young stellar objects (e.g., HK96, 
Sandell & Weintraub ISOUTI Abraham et al. l2U04] l. on the basis 
of statistical arguments it is believed that most low-mass pre- 
main sequence stars undergo FUor and/or EXor phases. Thus 
understanding the physics of the outbursts will shed light on 
fundamental processes of the early evolution of Sun-like stars 
and their circumstellar disks (structure, dust composition) and in 
turn on the changes of the conditions in the disks that govern the 
formation of planetary systems. The importance of the outburst 
phase in the early disk evolution was demonstrated for instance 
by Abraham et al. (|2009) , who discovered episodic crystalliza- 
tion of silicate grains on the disk surface due to the increased 
luminosity during the 2008 outburst of EX Lup, resulting in ma- 
terial that forms the building blocks of comets and planets. 

One of the best-documented and most-studied outbursts in 
the history of eruptive stars is that of V1647 Ori in 2003-2006 
(see the references collected by Aspin & Reipurth 2009 l. The 
low-mass central object (Mstar - 0.8 M©, Tstar - 3800 K, 
^star = 3.25 rJI Aspin et al., 1 20081 ), deeply embedded in the 
dark cloud LDN 1630, brightened by ~4.5 mag in Ic during three 
months (Briceno et al. 2004 ) until it reached a peak brightness in 
early 2004. It illuminated a new conical reflection nebulosity, the 
McNeil Nebula (McNeil l2004b . Subsequently, the object slowly 
faded during about 20 months. Between September-October 
2005 and early 2006 the object rapidly faded to its quiescent 
brightness (Acosta-Pulido et al. l20071 hereafter AP07). Variation 
of the near-infrared colors along the reddening path suggests 
that the brightening was partly caused by a temporal removal 
of an extinction slab (e.g., Reipurth & Aspin 2004). Optical 
and near-infrared spectroscopic observations (e.g., Reipurth & 
Aspin l2004l l showed characteristic spectral features attributed to 
accretion (Bry emission) and significant mass-loss (strong Ha 
emission with a PCyg profile). The variation of the accretion 
rate followed that of the optical brightness of the object, i.e., in- 
creased rapidly at the beginning of outburst, but decreased later 
on (AP07). The variation of the X-ray emission, which must 
originate from accretion processes, also follows the optical light 
curve, which supports the accretion scenario (Teets et al. 1201 11 1. 

Despite the numerous studies (see references in Aspin & 
Reipurth 2009 ), the circumstellar structure of V1647 Ori, and es- 
pecially its variations, have not been investigated in detail yej3- 
With the aim of studying the circumstellar disk and envelope 
of V16470ri and the dynamical processes during the eruption, 
we fitted the spectral energy distribution (SED) of the source 
at several epochs. In order to achieve a better characterization 
of VI 647 Ori, we initiated multi-epoch observations with the 
VLTI/MIDI mid-infrared interferometer. Spatial information on 



' The radius was recalculated from the post-outburst luminosity con- 
sidering a distance of d = 400 pc (Anthony-Twarog, T982)- 

^ The study of Muzerolle et al. ( I2005t was confined to using simple 
accretion-disk/remnant-envelope models for data of a single epoch. 



high angular scales was collected during two epochs. The first 
dataset from 2005 March was published in our preceding paper 
(Abraham et al. 120061 hereafter Paper I). Now we analyze these 
data together with the interferometric data from 2005 September, 
and with spectral energy distributions at different phases of the 
outburst. We performed detailed modeling of the circumstellar 
environment at each epoch to decide whether the observed tem- 
poral changes are related to the varying illumination of the disk 
by its central region or are caused by a change of the disk struc- 
ture. The results can be directly compared with those on another 
young eruptive star, PV Cep (Kun et al. 1201 11 1, whose observed 
flux variations are partly explained by variable extinction in the 
innermost part of the system, due to evaporation and reconden- 
sation of dust grains caused by the changing amounts of energy 
released during the outburst. Variations of the model parameters, 
together with their timescales, give insight into the dynamics of 
the outburst of the low-mass pre-main sequence star V1647 Ori. 

In Section 2 we describe the MIDI, Spitzer, and ground- 
based photometric observations and their data reduction. In 
Section 3 we present the results: the interferometric data, op- 
tical and infrared light curves, and compilations of the SED. In 
Sect. 4 we describe the radiative transfer models and discuss the 
temporal variation of physical parameters. Section 5 contains our 
conclusions. 



2. Observations and data reduction 

2.1. VLTI/MIDI 

During the 2003-2006 outburst, V16470ri was successfully ob- 
served twice with MIDI on the VLTI (Leinert et al. '2003'): on 
2005 March 2 and 2005 September 19. The projected baseline 
lengths were 56 and 62 m with position angles of 1 1 1 and 108°. 
These baselines allow one to marginally resolve structures of ob- 
jects with sizes of the mid-infrared-emitting regions of > 3 mas, 
which corresponds to > 1.2 AU at the distance of V1647 0ri. 
The two observing runs were carried out on the UT3-UT4 base- 
line of the VLTI. V16470ri was observed together with the 
spectrophotometric calibrator HD 37160. Owing to the optical 
faintness of the object and the lack of an adequate guide star, 
MACAO (the adaptive optics system of the VLTI for the Unit 
Telescopes) could not support the observations. At the time of 
our first observations, the seeing was w 1". During the observa- 
tions in September, it was better, 0.6 - 0.7". In this latter case, 
the atmospheric effects on the beam-overlap were less severe. 

The obtained sets of data consist of acquisition images with 
the N8.7 filter, 8 - 13//m low -resolution (R=30) spectra, and 
interferometric measurements. In the data reduction we fol- 
lowed the general processing scheme as described in previ- 
ous papers (e.g., Leinert et al. 2004 or Ratzka et al. 12007). 
MIDI data can be reduced in two independent ways: with the 
MIDI Interactive Analysis (MIA) package, which uses the power 
spectrum method, and the Expert Work Station (EWS) pack- 
age, which is based on a coherent, linear averaging method 
(Chesneau|2007B 

Since VI 647 Ori is not a bright object, we set 
up the observations in such a way that the scanning optical path 
delay (OPD) passes through the true zero OPD. In this case MIA 
should work more reliably than EWS. Nevertheless, we reduced 
the data with the 1.6 version of both softwares and found that 
the results of the two agree well. We also tested the results by 
using different mask widths, but found only a few percent effect 

^ The EWS-l-MIA package can be obtained from 
|http: //www. strw. leidenuniv.nl/~nevec/MIDI/ | 
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Table 1. Summary of observations of V1647 0ri. Notes give Spitzer AOR numbers and VLTI baseline parameters. Synthetic pho- 
tometry derived from the Spitzer spectroscopic data are marked with asterisks. The mid-infrared MIDI and Spitzer spectra are 
presented in Section 3. 



Date 


Instrument 


B and/W avelength 


Magnitude/flux [Jy] 


Notes 


References 


2007 Jan. 20 


IAC-80 


Ij 


19.24+0.06 mag 






2007 Oct. 28 


IAC-80 


Ij 


19.26±0.10mag 






2008 Feb. 6 


RCC 


Ic 


19.54 ±0.06 mag 






2007 Jan. 20 


TCS 


J 


14.60±0.06 mag 






2007 Oct. 22 


TCS 


J 


14.58+0.06 mag 






2007 Oct. 28 


TCS 


J 


14.44+0.05 mag 






2007 Jan. 20 


TCS 


H 


11.89+0.06 mag 






2007 Oct. 22 


TCS 


H 


11.96+0.12 mag 






2007 Oct. 28 


TCS 


H 


11.70+0.05 mag 






2007 Jan. 20 


TCS 


Ks 


10.01+0.04 mag 






2007 Oct. 22 


TCS 


Ks 


9.86+0.04 mag 






2007 Oct. 28 


TCS 


Ks 


9.77+0.06 mag 






2004 Mar. 28 


Spitzer/IRAC 


3.6 /jm 


1.883+0.058 


9463808 




2004 Oct. 27 


Spitzer/IRAC 


3.6 yum 


0.914+0.029 


12260864 




2005 Feb. 25 


Spitzer/IRAC 


3.6 yum 


0.945+0.034 


11570176 




2005 Mar. 25 


Spitzer/IRAC 


3.6 jim 


1.117+0.029 


11576320 




2004 Mar. 28 


Spitzer/IRAC 


4.5 yum 


2.864+0.087 


9463808 




2004 Oct. 27 


Spitzer/IRAC 


4.5 yum 


1.353+0.044 


12260864 




2005 Feb. 25 


Spitzer/IRAC 


4.5 yum 


1.443+0.051 


11570176 




2005 Mar. 25 


Spitzer/IRAC 


4.5 yum 


1.667+0.042 


11576320 




2004 Mar. 28 


Spitzer/IRAC 


5.8 yum 


3.713+0.111 


9463808 




2004 Oct. 27 


Spitzer/IRAC 


5.8 yum 


1.716+0.055 


12260864 




2005 Feb. 25 


Spitzer/IRAC 


5.8 yum 


1.839+0.062 


11570176 




2005 Mar. 25 


Spitzer/IRAC 


5.8 yum 


2.062+0.052 


11576320 




2004 Mar. 28 


Spitzer/IRAC 


8.0 yum 


5.721+0.172 


9463808 




2004 Oct. 21 


Spitzer/IRS (★) 


8.0 yum 


3.71+0.19 


12261120 




2004 Oct. 27 


Spitzer/IRAC 


8.0 yum 


2.745+0.090 


12260864 




2005 Feb. 25 


Spitzer/IRAC 


8.0 yum 


3.008+0.097 


11570176 




2005 Mar. 2 


VLTI/MIDI 


8.7 yum 


3.5+0.5 




Abraham et al. (2006) 


2005 Mar. 1 1 


Spitzer/IRS (*) 


8.0 yum 


2.32+0.12 


11569920 




2005 Mar. 24 


Spitzer/IRS (*) 


8.0 yum 


2.86+0.14 


12644096 




2005 Mar. 25 


Spitzer/IRAC 


8.0 


3.217+0.083 


11576320 




2005 Sep. 19 


VLTI/MIDI 


8.7 


2.0+0.5 






2004 Oct. 21 


Spitzer/IRS 


4 - 34 yum 




12261120 


Quanz et al. (2007) 


2005 Mar. 2 


VLTI/MIDI 


8-13 ytim 




56m/ 111° 


Abraham et al. (2006) 


2005 Mar. 1 1 


Spitzer/IRS 


4 - 34 yum 




11569920 


Quanz et al. (2007) 


2005 Mar. 24 


Spitzer/IRS 


4 - 34 yum 




12644096 


Quanz et al. (2007) 


2005 Sep. 19 


VLTI/MIDI 


8 - 13 yum 




62 m/ 108° 




2004 Mar. 15 


Spitzer/MIPS 


24 /jm 


18.2+2.8 


4320256 


MuzeroUe et al. (2004) 


2004 Oct. 14 


Spitzer/MIPS 


24 


16.5+2.5 


12260352 




2004 Oct. 21 


Spitzer/IRS (★) 


24 


13.68+0.71 


12261120 




2005 Mar. 1 1 


Spitzer/IRS (★) 


24 


8.65+0.45 


11569920 




2005 Mar. 24 


Spitzer/IRS (*) 


24 fim 


9.95+0.52 


12644096 




2004 Mar. 15 


Spitzer/MIPS 


70 fim 


25.7+1.7 


4320256 


Muzerolle et al. (2004) 


2004 Oct. 14 


Spitzer/MIPS 


lOfim 


30.7+1.9 


12260352 




2004 Oct. 15 


Spitzer/MIPS (★) 


70 fim 


36.6+3.7 


12260608 




2005 Mar. 1 


Spitzer/MIPS (★) 


70 fim 


30.6+3.1 


11570432 




2005 Mar. 4 


Spitzer/MIPS (★) 


70 fim 


32.6+3.3 


11576576 




2006 Nov. 8 


Spitzer/MIPS 


70 fim 


2.7+0.2 


17455360 




2004 Oct. 15 


Spitzer/MIPS 


55 - 95 yum 




12260608 




2005 Mar. 1 


Spitzer/MIPS 


55 - 95 yum 




11570432 




2005 Mar. 4 


Spitzer/MIPS 


55 - 95 yum 




11576576 





on the visibilities. Although we investigated all calibrator data 
taken on the nights of our MIDI observations, one should proba- 
bly not consider those where MACAO was used to estimate the 
errors of the instrumental visibilities. Note that MACAO was not 
used for the observations of HD 37160. Nonetheless, the results 
for HD 37160 agree well with the other instrumental visibilities 
(within Icr) obtained the same night. The beam overlap was good 
at both epochs, even despite the lack of MACAO support. We 
found that the eiTors of the correlated fluxes are small (a few 



percent), but those of the total spectra are higher (10-15%) and 
can change slightly with data reduction parameters. For the sake 
of simplicity, we applied a 15% error on the spectra, independent 
of wavelengths. The errors of 15% on visibility data (see Fig.|2]i 
are consistent with other observations (c.f. Chesneau l2007l l. 

Similarly to Paper I, we derived 8.7/im photometry of 
VI 647 Ori from the acquisition images, using the acquisition im- 
ages of HD37160 for calibration (see Table[TJ. 
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In Paper I we presented the first dataset from 2005 March 
reduced with MIA 1.3. For the sake of consistency, we re- 
reduced the dataset with MIA 1.6. We discuss the differences 
between the results produced by the different versions of MIA in 
AppendixlAl 

2.2. Optical and near-infrared pliotometry 

In addition to the already existing photometric data (AP07), new 
near-infrared (NIR) J-, H-, and Ks-band observations were car- 
ried out in 2007 and 2008 (Table [T]i. These data can be used 
to characterize the post-outburst quiescent phase object. The 
data were obtained using CAIN-2 installed on the 1.52 m Carlos 
Sanchez Telescope (TCS) at the Teide Observatory (Tenerife, 
Canary Islands, Spain). We also observed V1647 0ri in the Ij 
band with the 82 cm IAC-80 telescope at the Teide Observatory 
and in the Ic band with the 1 m RCC telescope at the Piszkesteto 
station of the Konkoly Observatory. The technique of observa- 
tion, data reduction, and photometric calibration of these data 
were identical with our previous TCS, IAC-80, and RCC obser- 
vations of V1647 Ori and are described in detail in AP07. 

2.3. Spitzer arcliival data 

VI 647 Ori was observed with the Infrared Array Camera 
(IRAC) onboard the Spitzer Space Telescope at 3.6, 4.5, 5.8, and 
8.0 yum at four different epochs (Table[T]i using the subarray read- 
out mode. Offering exposure times as short as 0.02 s, the subar- 
ray mode enables the observation of bright sources that would 
saturate the detector in other readout modes. IRAC data can be 
downloaded from the Spitzer Heritage Archive. We considered 
all IRAC sub array images that provided more reliable results 
than the mapping data. These were obtained on 2004 March 7 
and the results were published by Muzerolle et al. (2(H)5 l. The 
authors derived somewhat 10%) lower fluxes at each observ- 
ing wavelength. All these observations except for those carried 
out in 2004 October consist of nine different dither positions 
with 64 images at each dither position (with size of 32x32 pixel). 
They were obtained with an exposure time of 0. 1 s per frame. 
Four dither positions were used in 2004 October. The data pro- 
cessing started with the Spitzer Science Center (SSC) basic cal- 
ibrated data (BCD) produced by the pipeline version SI 8.7. The 
BCD image cubes of the 64 frames were combined into two- 
dimensional images using the irac-subcube-collapse IDL routine 
provided by the SSC. Aperture photometry was performed on 
the final images at each wavelength using a modified version of 
the IDLPHOT routines. The aperture radius was set to 3 pixels 
(pixel scale 1'.'2), the sky background was computed in an an- 
nulus with an inner radius of 3 pixels and a width of 4 pixels. 
In the course of the sky estimates we used an iterative sigma- 
clipping method, where the clipping threshold was set to 3cr. 
Following the outline of Hora et al. (2008), we also applied an 
array-dependent photometric correction and a pixel-phase cor- 
rection to the measured flux densities. An aperture correction 
was then performed using the values published in the IRAC in- 
strument handboolfl (IIH). The flux density values measured at 
the different dither positions (at each band) were averaged to ob- 
tain the final photometry. The final uncertainties were computed 
by quadratically adding the measurement errors (obtained from 
the individual flux density values in each band) and an absolute 
calibration error of 3% (IIH). No color correction was applied 



since it would have only a marginal effect (< 1%) according to 
our calculations. 

VI 647 Ori was observed with the Infrared Spectrograph 
(IRS) of Spitzer at three different epochs. On 2004 October 21 
and on 2005 March 24 the target was measured with the short- 
low (5.2 - 14.5 yum, R~60 - 127), short-high (9.9 - 19.5 //m, 
R~600) and long-high (18.7 - 37.2 //m, R~600) channels, 
while on 2005 March 1 1 the short-low and the long-low (14.0 - 
38.0 fim, R ~ 57 - 128) channels were used. All these data were 
published by Quanz et al. (I2007I I. and we considered these spec- 
tra for our modeling. 

VI 647 Ori was also observed with the Multiband Imaging 
Photometer for Spitzer (MIPS, Rieke et al. 2004) in scan-map 
mode on 2004 March 15 and in photometry-mode on 2004 
October 15 and 2006 November 8. The data analysis of the MIPS 
observations started with the pipeline (S16.1) generated basic 
calibrated data (BCD) images. At 24 /vm we performed a flat- 
field correction and a background-matching using the MOPEX 
tool (Makovoz & Marleau 2005). Following the general steps for 
processing MIPS 70yum data described by Gordon et al. (2007), 
we performed column-spatial-filtering and time-median-filtering 
on BCD images obtained at this wavelength. The improved BCD 
images were subsequently mosaiced using MOPEX. The output 
mosaic was resampled to 2'.'45 and 4"pixel"' at 24 and 70jum. 
The MIPS 70 fim array is affected by nonlinearity at high count 
rates. To correct for the 70 //m nonlinearity effects, we applied 
the formula described by Dale et al. (120091 see Eq. 1). At 24 fim 
VI 647 Ori is saturated, therefore we used a model PSF to fit 
the pixels that were still in the linear regime. At 70 fj.m we used 
a modified version of the IDLPHOT routines to obtain aper- 
ture photometry for VI 647 Ori. The aperture was placed at the 
SIMBAD position of the source and the background was esti- 
mated in a sky annulus between 39"and 65". For images taken at 
the first two epochs the aperture radius was set to 18". In the last 
image, however, the source appeared at the edge of the image, 
which allowed us to use only a smaller aperture (radius of 12"). 
The aperture-correction factors, appropriate for a 60K black- 
body, were derived following the outline of Gordon et al. (2007). 
The final uncertainties were computed by quadratically adding 
the internal error and the absolute calibration uncertainty (we 
adopted 7%, see MIPS data HandboolQ). The final photometry 
and its uncertainty is listed in Table [1] The MIPS data in Table [1] 
are color-corrected (by a factor of 0.91). Our results for 2004 
March are higher by 16% and 46% at 24 and 70 fim, compared to 
those of Muzerolle et al. (I2005I I because of the different pipeline 
versions and the additional data reduction steps. 

VI 647 Ori was also observed using the MIPS SED mode. 
Low-resolution far-infrared (55-95jum; A/AA ~ 15-25 ) spectra 
were obtained on 2004 October 15 and on 2005 March 1 and 4. 
The data reduction of SED observations also started with the 
BCD images (pipeline version S16.1) and MOPEX was used for 
the processing steps (combination of data, background removal, 
application of the dispersion solution) and to compile the final 
images with a pixel scale of 9.8". An aperture covering 5 pixels 
was used to extract the spectra from the final images. The aper- 
ture correction factors were taken from Lu et al. (12008) . 

Finally, synthetic 8.0, 24, and 70 fim photometry was derived 
from the IRS and MIPS SED data by convolving the spectra with 
the corresponding IRAC or MIPS filter profiles. The results are 
presented in Table [T] 



Version 1.0, http://ssc.spitzer.caltech.edu/irac/ 
iracinstruinenthandbook/ 



' Version 2.0, http://ssc.spitzer.caltech.edu/mips/ 
mipsinstrumenthandbook/ 
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Fig. 1. Light curves of VI 647 Ori. The vertical dashed Unes de- 
note the epochs of the MIDI observations. Iq (O.Sjum) and Ks 
(2.2//m) band data are from AP07 and this paper. The Ic data 
are scaled up by a factor of 40 for the sake of better plotting. 
S.Ojum data are from MuzeroUe et al. (2005 1, Paper I, Quanz et 
al. 12007), this paper (obtained from MIDI acqusition data), and 
Aspin et al. (2008, extrapolated from their N' band data). 24 fim 
data are from this paper and Aspin et al. ,2008. (extrapolated from 
their band data). 70 fj.m data are from this paper. 



3. Results 

3.1. Optical and infrared liglit curves 

Figure [1] shows the light curves of VI 647 Ori at five optical 
and infrared wavelengths (0.8, 2.2, 8.0, 24, and IQfim) be- 
tween 2004 and 2007. The overall shape of the Ic light curve, 
which is the most complete, can be divided into a plateau (2004 
February - 2005 September), a rapid fading (2005 October - 
2006 February), and the subsequent quiescent phase. In addi- 
tion to the general slow fading during the plateau phase, short 
timescale low-amplitude variations are also observed. The rapid 
fading started around the epoch of the second MIDI observa- 
tions. The available data in the Ks band may suggest a light curve 
similar to that of Ic, but with a smaller amplitude. Although the 
light curves at 8, 24, and 70yum are even more sparsely sam- 
pled, they are not inconsistent with the shape of the Ic band light 
curve. We note that at 70 //m, the object became ~ 20% brighter 
between the 2004 March peak and 2004 October. During the 
same period, the IRAC fluxes decreased by a factor of 2. This 
suggests a time shift between the peak brightnesses at the differ- 
ent wavelengths. The cause of such a shift at 70//m is unclear. 

Comparing the flux values obtained in 2006 to the pre- 
outburst ones (interpolated in wavelength in Fig. 1 of Abraham 
et al. 2004), we conclude that V1647 Ori returned to quiescence 
after 2005. The 2007 fluxes at 8 and 24 //m are even lower than 
the respective pre-outburst values. This could either be due to 
variability in the quiescent phase (similar as seen in the pre- 
outburst Ic data of Briceno et al. 2004), might be related to the 
beam differences of the different instruments, differences in the 
photometric systems, or might even indicate a different object 
structure after the outburst. 



3.2. Mid-infrared interferometry 

Fig. |2] shows the calibrated visibilities obtained in 2005 March 
and September. The baselines of the two observations are very 
similar (Table [T]i, thus the two visibility curves can be directly 
compared. The difference indicates a significant change in the 
geometry of the mid-infrared emitting region of VI 647 Ori in 
the six months that elapsed between the two observations. As 
a first approximation, we considered Gaussian brightness dis- 
tributions to estimate the characteristic size of the mid-infrared 
emitting region defined as the FWHM of the Gaussian. Although 
the object is much more resolved at the second epoch, and ac- 
cordingly the Gaussian assumption is less adequate for size esti- 
mation, we converted the visibilities to the FWHM, which pro- 
vided a simple comparison of the interferometric data. Single 
Gaussians with a FWHM=10 and 16 mas, equivalent to 4.0 and 
6.4 AU at the distance of 400 pc, respectively, fit both visibil- 
ity datasets well longward of ~ lOjum. By fitting Gaussians to 
the visibilities, characteristic sizes of 5-7 mas (2.0-2.8 AU) and 
1 1-13 mas (4.4-5.2 AU) can be estimated at shorter wavelengths 
for the two epochs. The object was less resolved at wavelenghts 
shortward of ~ 10/vm at both epochs, probably due to an in- 
ner warm part of the system. Because the shape of the visibility 
curves is similar to those observed for other young stellar objects 
(YSOs) with MIDI (e.g., Leinert et al. 2004, Quanz et al. i20()6l 
Ratzka et al. 2007), we can assume that they have similar cir- 
cumstellar structures, which means that the radiation from the 
inner edge of the disk can dominate at short wavelengths. 

Similarly to Paper I, we do not see any sinusoidal variations 
in the spectrally resolved visibilities that are potentially caused 
by a companion. Since the baselines of the two observing runs 
are very similar, we can only repeat our earlier finding that the 
shape of the new visibility curve suggests that no companion 
whose separation is between about 10 and a few hundred AU and 
has a brightness ratio higher than 10% is present at the measured 
position angle . 

Our two 8 - 13//m MIDI spectra are shown in Fig|2] 
Although the absolute level of the mid-infrared spectra de- 
creased by a factor of 2 between the two epochs, the shapes 
of the two spectra are similar. In Fig. |2] (d) we plot the first 
MIDI spectrum together with the two IRS spectra obtained in 
2005 March. Variations on weekly timescales are clearly present. 
Quanz et al. (2007) found weak silicate emission in the IRS spec- 
tra. Although the ozone band can affect the spectra at 10/im, this 
weak feature might also be present in the MIDI spectra. 

The uncorrected spectra are calculated as the difference of 
the total and the correlated ones. The correlated and uncorrelated 
spectra of an object are dominated by radiation from different 
circumstellar regions: the inner regions of a few AU size and the 
outer parts (van Boekel et al. I2004I I. In the case of V1647 0ri, 
most of the mid-infrared emission originates from the inner com- 
pact zone of the circumstellar environment, the correlated spec- 
trum being higher than the uncorrelated spectrum. The corre- 
lated fluxes decreased significantly, by approximately 60% be- 
tween the two MIDI observations, while the relatively low un- 
correlated fluxes increased by about 30% longward of 10/im, 
and by up to 70% shortward of 10 fim (Fig.[2](f)). These findings 
may suggest either a structural change in the circumstellar envi- 
ronment of the object or that the fading of the inner and the outer 
parts of the mid-infrared emitting regions of VI 647 Ori occured 
on different timescales. The correlated spectra look featureless; 
the weak silicate emission seems to be associated with the outer 
parts of the system. 
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Fig. 2. Calibrated visibilities (a) and Gaussian size estimates (b). MIDI spectra are shown in panel (c). The MIDI and Spitzer/IRS 
spectra from 2005 March are compared in panel (d). The errors of the Spitzer spectra (~10%) are not shown for the sake of clarity. 
Correlated spectra (i.e., inner disk spectra, panel (e)) and uncorrected spectra (i.e., outer disk spectra, panel (f)) of V1647 Ori are 
also shown at two epochs. The baselines of the two observations are very similar, so the MIDI data of the two epochs can be 
compared. The shaded areas mark the errors in the panels (a,c,e,f). The vertical dotted lines show the ozone band that could affect 
the ground-based observations. 



3.3. Spectral energy distribution 

We constructed SEDs of V1647 0ri for several different epochs 
to study its circumstellar environment. These epochs were 2004 
March (peak of the outburst), 2004 October (plateau phase 1), 
2005 March (plateau phase 2), 2005 September (shortly before 
the rapid fading). The latter two epochs are those of the MIDI 
observations. Due to the lack of simultaneous optical and near- 
infrared data, we included data points in the SEDs obtained 
within a period of a few weeks of the nominal date. We also com- 
piled a combined SED from pre- and post-outburst quiescence 
data (cf. Aspin et al. l2008l l. Because of the sparse time coverage 
of observations, i.e., no simultaneous SED, we merged all avail- 
able measurments and plotted all pre- and post-outburst data to- 
gether in the quiescent SED. The five SEDs are plotted in Fig.|3] 
The Spitzer IRS and MIPS SED spectra are also plotted. The 
absolute errors of these spectra are ~ 10%. For reasons related 
to the radiative transfer modeling, which is described in Sect. 21 
we consider data only shortward of 100 fim. Submillimeter and 
millimeter data are used only for the estimation of the mass of 
the system (see Sect. 4.1.1). 

4. Modeling 

In the following we fit both the SEDs and the visibilities using 
a radiative transfer code to provide a model of the circumstellar 
environment and the heating mechanisms. Our strategy for fit- 
ting the data was as follows. First we set up a reference model 
that fits the 2005 March data (Sect. We chose this epoch 



to serve as a starting point because our observational data set is 
most complete here. Then we modified this initial setup as neces- 
sary to fit the data of the other epochs (Sect. l43t . The significant 
change of the MIDI visibilities might be explained by different 
scenarios, which are discussed separately in Sect. 14.41 



We note that first we tried to adopt the simplistic disk model 
of the VI 647 Ori system presented in our Paper I for all epochs. 
Although the SED models were satisfactory, the model cannot 
reproduce the low visibility values measured in September 2005. 



4.1. General overview of the model 

For the radiative transfer modeling we used the Monte Carlo 
code MC3D (Wolf et al. 1999. Wolf l2003l Schegereret al. i2008l l. 
We assumed an axial symmetry of the circumstellar environ- 
ment, therefore we used a two-dimensional geometry in polar 
coordinates (r, 9). 

The model system consists of a central star, an inner accre- 
tion (gas) disk, and a dusty outer disk embedded in an envelope, 
from which a cone-shaped cavity is cut (Fig.lUi. For the structure 
of the disk we assumed the density profile to be like the one used 
by Wood et al. (,2002i i. 
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Fig. 3. Spectral energy distributions of V16470ri at five different epochs. Because of the short timescale variability at the optical 
and near-infared wavelengths, we present photometry from between 2004-02-01 and 2004-04-30 (2004 March), 2004-10-01 and 
2004-10-31 (2004 October), 2005-02-01 and 2005-04-30 (2005 March) and 2005-09-03 and 2005-10-03 (2005 September) in these 
figures. Our models to the SEDs, where the inner radius of the dust disk is set to 0.5 AU in quiescence and 0.7 AU in all outburst 
epochs, are overplotted with solid lines (see Sectj4). Between the outburst epochs only the accretion rate was changed here. Note 
that the model for 2005 September shown here does not fit the MIDI data (see Fig.|5]and Sect. 14.41 ). The last panel shows all models. 
The decrease of the brightness of the object with time, from peak (top) to quiescence (bottom) is continuous. References: quiescence 
- Abraham et al. (2004), Paper I, Aspin et al. (2008i; peak (200 4 Mar ch) - AP07, McGehee et al. (2004), Muzerolle et al. (2005); 
plateau phase I (2004 October) - AP07, this work, Quanz et al. (120071 ): plateau phase II (2005 March) - AP07, this work, Quanz et 
al. ( 120071) : before the rapid fading (2005 September) - Aspin & Reipurth (120091 ), this work. 




gas dusty disk disk 
disk & envelope 



Fig. 4. Schematic picture of the model geometry of the object 
(not to scale). The line of sight crosses through the envelope. 
The inclination of the disk is i? 60°, the opening angle of the 
conical cavity is 50° (AP07). 



where m is the radial coordinate in the disk midplane, z is the 
vertical height and po is determined from the mass of the disk. 
The scale height him) increases with radius as, 

''^'"(lo^f' 



where ho is the scale height at a radial distance of 100 AU. The 
density profile of the envelope is spherically symmetric: 

Pem ~ PO.env 1 j 

Here r gives the distance from the star, and where po.env is de- 
fined as /p(),disk at a certain arbitrary point of the disk and the 
envelope. The value of the constant / has no physical mean- 
ing, it is set so as to separately control the mass of the disk and 
the envelope. The transition between these two components is 
assumed to be smooth. In vertical direction the envelope starts 
where the density of the disk decreases to the density of the en- 
velope. In the cavity we used a density distribution identical to 
that of the envelope, but less dense by a factor of 10"^. 

The temperature distribution of the system is determined by 
the heating sources: the central star, the heated dust grains, emit- 
ting blackbody radiation, and the accretion, consisting of an a- 
type disk concentrated on the midplane and a hot spot on the 
stellar suface (for details of accretion effects see Schegerer et 
al. i2008t . The accretion process is characterized by three param- 
eters, the accretion rate M, the temperature Tspot of the accreting 
region on the surface of the star, and the magnetic truncation 
radius Rtmnc- After calculating the temperature distribution, the 
SED and a projected image of the system is produced at an in- 
clination angle § with a ray-tracer. 

The circumstellar disk in our model is embedded in a signif- 
icantly larger envelope that stretches out far beyond the disk, to 
3000 AU, in accordance with the mm-images from Tsukagoshi 
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et al. (I2005I I. The inclination of the disk was fixed tohe-& x 60° 
and we set the opening angle of the cavity to 50° (AP07), thus the 
line of sight to the star crosses the envelope. Although the mid- 
infrared ice features indicate the presence of foreground material 
(Vacca et al. '2004' Rettig et al. 120051 Quanz et al. l2007b . the ice 
can also be in the outer regions of the extended envelope around 
VI 647 Ori. Therefore we fixed the interstellar extinction to be 

= and assumed that all extinction is circumstellar. 

The lOjum silicate feature is too weak to determine the dust 
composition via model-fitting. Therefore we considered stan- 
dard interstellar dust, a mixture of 62.5% astronomical silicate, 
and 37.5% graphite (the optical parameters and grain size distri- 
bution were taken from Weingartner & Draine 2 001) . The shape 
of the continuum-subtracted features (Quanz et al. I2007I I and the 
presence of an envelope support the plausibility of this assump- 
tion. We used spherical dust grains, with a power-law size dis- 
tribution n(a) oc a^^-^, for flmin < fl < Qmax, where n{a) is the 
number of dust particles with a radius a. For the minimum and 
maximum grain size we used 0.005 ;um and 0.25 fim. To limit the 
number of free parameters in the model fitting, we fixed a^in and 
Qmax- Note that increasing the maximum grain size to a few yum 
would decrease the silicate emission, but would not change the 
model SEDs significantly. The gas-to-dust mass ratio is assumed 
to be 100, while the grain mass density is set to pg -2.1 gcm"^. 
We have found that with the selected grain size distribution we 
underestimate the submm-to-mm flux in all our models, unless 
we assume a disk mass of » IM©, which would result in gravi- 
tational instabilities. In addition, the submm - mm slope of the 
model is steeper than that of the measured data, which leads to 
the conclusion that the system possibly contains larger grains 
than the maximum size used. However, we think the larger grains 
are within the inner regions of the disk closer to the midplane 
due to sedimentation, affecting the SED significantly only in the 
(sub-)mm regime, since this region is optically thick at shorter 
wavelengths. Consequently, we consider our model to be valid 
shortwardof ~ 100/zm. 

We estimated the mass of the system from the 1.3 mm 
flux (Lis et al. 11999). using the method described in Beckwith 
et al. ( fT999] l. Assuming T = 50K and = 0.02 cm^g ', 
we found that the total mass of the circumstellar matter is 
Mtot - 0.045 Mq. The estimate is similar to that of Andrews 
et al. ( 120041) . We considered this estimate as a fixed parameter 
during our model-fitting. 

4.2. The reference model: 2005 March 

As a first step we fitted the available photometric and interfer- 
ometric data for the 2005 March epoch. Since the SED is not 
composed of results of simultaneous observations and the ob- 
ject is variable on short timescales (see Sect. 13.31 1. we consider 
ranges of data at certain wavelength^ Such a dataset prevents 
us from calculating correct meaningful values and compar- 
ing diff'erent models in a mathematical way. Instead, we eval- 
uated the plausibility of the models by visual inspection of 
the fits - as done in most of the cases, see e.g., Schegerer et 
al. ("Sods'), Ratzka et al. (l2009b . van Boekel et al. (l20T0l) . Juhasz 
et al. (2012 ). The model-fitting in the latter two cases is very sim- 
ilar to our case because the objects studied by the authors (T Tau 
and EXLup, respectively) are also variable and the authors fitted 
merged SED datasets, as we did for the different epochs. 

We set the derived model parameters against results found 
in the literature to filter out unphysical models that still showed 



good fits. The model has many parameters, and a good fraction 
of these cannot be determined or constrained by measurements. 
We attempted to fit the data with the extreme cases of envelope- 
only and disk-only models, but we conclude that neither of the 
two components can be neglected. The parameters of the cir- 
cumstellar environment are not clearly defined by the models 
because the parameters of the envelope and the disk are not in- 
dependent (e.g., Natta 1993 ). To reduce the degeneracy we fixed 
additional system parameters and tried to consider canonical val- 
ues. We set the outer radius of the disk to 500 AU and the value 
of Q-env to -1.5 (Shu |1977] l. To reduce the number of free pa- 
rameters even more, we set the inner radii of the disk and the 
envelope to be identical. The investigated parameter ranges are 
listed in Table|2l In the following we present our results with the 
fixed parameters, although our conclusions and the trends seen 
- based on the models we propose - are also valid for slightly 
different setups of the envelope-disk system. Therefore we do 
not present best-fit models here, but a series of plausible good 
models. 

During the fitting procedure different M and Rin values were 
considered, with step sizes of 10"^ Mo/yr and 0.05 AU, respec- 
tively. The step sizes were 5 - 10% maximum for all other 
epochs. For the reference epoch we found that the accretion rate 
is 3.5 X 10"^ Mo/yr and the inner radii of the disk and the enve- 
lope are R^^ - 0.7 AU. Since the line-of-sight to the star crosses 
the envelope but not the disk (Fig. 4), the parameters of the en- 
velope (mass, Rin and aenv) strongly determine the fits in the 
optical-NIR regime. The thick envelope increases the extinction 
in the system in order to decrease the otherwise overestimated 
optical-NIR model fluxes to the measured valuefl and it also 
helps reducing the height of the intrinsically strong 10 /xm sili- 
cate emission feature. The parameters of our model are given in 
Table|2] the model can be seen in the upper panels of Fig.|5] The 
visibilities are somewhat underestimated. A better fit could be 
achieved, but at the cost of a poorer SED-fit. 

Comparing the disk and envelope parameters of our model 
with the results of similar modeling efforts of embedded YSOs 
in the literature (e.g., Whitney et al.'2003), we may conclude that 
the circumstellar environment of V1647 Ori is rather typical, and 
no structural features directly linked to the eruptive nature of the 
object can be identified. With our model we could describe the 
system with a parameter setup that is plausible for young stellar 
objects. 

4.3. Modeling the variations during outburst 

In order to fit the SEDs of the other epochs, we used the 2005 
March model as an initial setup. At first sight, the SEDs differ 
only in their absolute levels, but their shapes are very similar 
(Fig. |3]). Therefore, we attempt to adjust the shift between the 
levels by the variation of the accretion rate. This way we test if 
it is possible to keep the geometry of the circumstellar environ- 
ment identical, and only change the illumination of the system 
as the outburst proceeds. Additionally, we tried to vary only the 
smallest number of parameters possible between the different 
epochs and restricted ourselves to applying modifications that 
could be the result of a realistic process and may occur in such 
a short period of time. Significant changes, affecting a large part 
of the disk or a complete restructuring, could happen only on 
longer timescales (Chiang & Goldreich |1997l l. 

According to the standard picture of eruptive YSOs (HK96), 
the accretion is the main energy source in the system in outburst. 



* This holds also for the other epochs. 



' Note the absence of NIR photometric data in 2005 (Figs.[T]and[3j. 
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Table 2. Parameters of the best-fitting model for 2005 March. The fitted parameters are shown in italics. 
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Parameters 


Final model 


References for fixed parameters 


Ranges for fit 


Stellar parameters 








Temperature (T^f^y) 


3800 K 


Aspin et al. <20Q8|l 




Mass (Mstar) 


0.8 Mo 


Aspin etal. (20081 




Radius (R,,.^) 


3.25 Ro 


Asoin et al. (i2008il. recalculated 




Interstellar visual extinction (Ay) 


mag 






Circumstellar disk parameters 








Inner radius of dusty disk (Rm,Asv) 


0.7 AU 




0.2- 1.5 AU 


Outer radius of dusty disk (i^out.disk) 


500 AU 






Scale height at 100 AU (Ho) 


15 AU 




3 - 25 AU 


Flaring index (fi) 


1.2 




I-'/, 


Exponent of radial density profile (ffjisk) 


-1.75 




-2.25 - -1.5 


Total mass of disk and envelope (M) 


0.045 Mo 


estimated 




Distance (d) 


400 pc 


Anthony-Twarog (|1982|, AP07 




Inclination (ff) 


60° 


AP07 




Circumstellar envelope parameters 








Inner radius of dusty envelope (-Rin.ciw) 


0.7 AU 




0.05- 1.0 AU 


Outer radius of dusty envelope (-Rout.cnv) 


3000 AU 


Tsukaeoshi et al. (120051 




Exponent of radial density profile (ttenv) 


-1.5 






Parameters for the accretion 








Accretion rate (M) 


3.5xI0-«Moyr-i 




2- lOxlO^^Moyr-i 


Magnetic truncation radius (-/?tmnc) 




Calvet Sl Gullbring (|l998l 




Temperature of the hot spot (Tspot) 


6500 K 


Calvet & Gullbrine Jl998ll 





while its contribution in quiescence is less pronounced. The ac- 
cretion rate is expected to increase from the pre-outburst phase 
to the peak of the outburst and decrease during the eruption as 
the source is fading. The modeling results confirmed this ex- 
pectation. We found that the accretion was strongest at the peak 
brightness of 2004 Mai'ch (7 x 10"^ Mg/yr, see Tab.©, and grad- 
ually weakened as the source approached the end of its flare-up, 
but still accreted strongly in quiescence (3 x 10"^ M^lyr). The fit 
to the mid-infrared part of the SED depends mainly on the accre- 
tion rate. The change found in the accretion rate is solid, since 
the fading can be seen clearly at these wavelengths (see Fig.|3]l. 
The strategy of changing the accretion rate worked for the whole 
SEDs of the outburst epochs, but failed to reproduce the second 
MIDI interferometric data from 2005 September (see Sect. 14.41 
and the middle lower panel of Fig. |5]l. Lowering the accretion 
rate to ^ 5% of the peak value reduced the fluxes considerably 
at all optical-IR wavelengths in quiescence. This model fits the 
mid-infrared part of the SED well, but the excess of the model in 
the optical-NIR part of the SED indicates that the accretion rate 
is not the only parameter that must be adjusted. 

The inner radii of the dusty components in the reference 
model are set to 0.7 AU, which provides good fit for the other 
two outburst epochs in 2004. Accordingly, the geometry of the 
inner part of the system was unchanged during the plateau phase. 
For quiescence, we found the best fit by using /?in,qui = 0.5 AU, 
which is the radius of both the dust disk and envelope. The vari- 
ation of the inner radii of the disk and the envelope indicates 
that a dynamical process worked immediately after the outburst, 
which enlarged the inner dust-free hole from 0.5 to 0.7 AU on 
a timescale of a few months. The derived = 0.5 and 0.7 AU 
radii correspond to lOOOK and 1500 K in quiescent state and 
at the peak brightness phase, respectively. These temperature 
values are close to the canonical dust sublimation temperatures 
(~ 1300 - 1600 K). However, moving the inner radii of the qui- 
escent model inward to the corresponding sublimation radii, the 
model underestimates the optical data points due to the increased 
extinction. Although similar deviations of the radii were found 
at other eruptive young systems (e.g., Sipos et al. 120091 Juhasz 
et al. 120 12i Eisner et al. 1201 11 1, the cause of this larger inner ra- 



dius is not clear. Varying parameters of the models are shown in 
Table [3l 

4.4. 2005 September: The puzzling MIDI data 

For the 2005 March and September epochs the MIDI visibili- 
ties give additional constraints for the system parameters. We at- 
tempted to fit the 2005 September SED and visibility values us- 
ing the same method followed above. We kept the inner radii un- 
changed (Rin = 0.7 AU ) and changed the accretion rate between 
the two MIDI epochs from 3.5 x 10 MJyi to 2.0 x 10"^ Mo/yr. 
Due to the fading of the central source the mid-infrared emit- 
ting region is smaller and results in slightly increased visibilities. 
This is just the opposite of what the measurements show (shown 
with open triangles in Fig. |5] lower middle panel). Clearly, al- 
though for the SEDs our first strategy of changing only the ac- 
cretion rate and the inner radii of the dust disk and the envelope 
works, it cannot account for the changes in the visibilities. 

In the framework of our static radiative transfer models the 
geometry of the circumstellar environment also had to undergo 
changes in the inner regions of the disk-envelope system to ex- 
plain the lowering of the visibility points. The change has to take 
place in the inner regions of the system, because the dynamical 
timescales for changing the structure of the environment farther 
out is much longer than the half-year elapsed between the two 
MIDI observations. 

What that might help to explain the changes in the visibility 
data is to assume that due to the increased heat - because of the 
enhanced accretion - the scale height of the regions close to the 
inner edge of the disk increases in the first phase of the outburst. 
By introducing such a region to the system and varying the scale 
height of the inner part of the disk, i.e., a variable pujfed-up in- 
ner rim, we also see variations in the visibilities. The visibilities 
are higher if the scale height increases. A higher inner rim means 
an additional hot component in the system, which increases the 
visibilities. However, this rim also casts more shadow on the 
disk, thus certain regions behind the rim will not be illuminated, 
which might also affect the visibilities. We tested several geome- 
tries of the disk to see how strong the effect is and found that this 
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Fig. 5. Our models for the 2005 March and September epochs are shown in the top and bottom row. The SEDs and the calibrated 
visibilities are shown in the left and middle columns. Model visibilities were calculated at eight different wavelengths. In the bottom 
left panel the solid line, in the bottom middle panel the open triangles show the model where only the accretion rate was decreased 
relative to the reference model. This model fits the SED quite well, but not the visibilities. The dash-dotted line {bottom left panel) 
and the filled triangles {bottom middle panel) show the model where the inner radius of the envelope was also increased. As an 
alternative interpretation of the mid-infrared fits, the observed mid-infrared total, correlated and uncorrelated MIDI spectra (dashed, 
solid and dotted lines, respectively) and the corresponding model values (filled circles, diamonds and downward pointing open 
triangles, respectively) are shown in the right column. The lower right panel refers to the case where the inner radius of the envelope 
increased well beyond the dust sublimation radius, inducing larger changes in visibilities than the ones described earlier. /?in,disk was 
kept unchanged. Data probably affected by the atmospheric ozone layer are not shown. 



process could only account for approximately one quarter of the 
observed decrease, i.e., 0.05-0.08 in the visibilities. Therefore 
we do not discuss these models further. 

Another straightforward idea is to increase the inner radii of 
the dusty components from March to September, which we ex- 
pect to have the effect of increasing the mid-infrared emitting 
region, i.e., making the source more resolved. The inner radii of 
the envelope and the disk were changed previously and seem to 
be the parameters that could vary the easiest (compared to disk 
scale height or density profile exponents). However, any trans- 
portation mechanism that removes material is much more effec- 
tive on the significantly less dense envelope, and would move its 
inner radius much farther out. Therefore we do not consider any 
change in the disk in the following. 

In our model that fits the visibilities, the dusty envelope starts 
only at 3.0 AU, and the matter removed from the inner regions 
is piled up in a thin shell close to the inner edge of the envelope. 
This distance corresponds to ^ 600 K temperature. By moving 
the inner edge outwards, the infrared flux increased slightly, so 
we had to adjust the accretion rate again, i.e., lowered to 1.6 
Moyr"' X 10 Moving matter outwards, and creating a shell of 
matter has the additional effect of reducing the extinction, thus 
increasing the optical-near infrared flux. As the results show in 
Fig. |5] improving the fit to the visibilities deteriorated the SED- 
fit. The extinction in this model dropped to Ay =11.5 mag, which 
would have led to an unobserved brightening of VI 647 Ori (cf. 
Fig. 2). Although the accretion luminosity decreased in this pe- 
riod, it could hardly compensate for the effect of AAy. The un- 
certainty of the optical part of our model for 2005 September 



prevents us from analyzing of the evolution of the extinction at 
this epoch. We attempt to find the mechanism that can produce 
the dramatic variation of the envelope radius between March and 
September 2005 in Sect. Ill 

5. Discussion 

5.1. The varying accretion rate 

The peak accretion rate value we found is on the same order of 
magnitude as the 5 - 10 x 10"^ Mg/yr values derived from Bry 
measurements (AP07) and modeling (Muzerolle et al. 2005). 
During the plateau phase of the outburst, as the object slowly 
faded, our model accretion rate slightly decreased. Both obser- 
vations and modeling indicate that x 10"^ M© material was ac- 
creted onto the star during the outburst. 

The derived accretion rates fit in the canonical picture of 
eruptive low-mass YSOs. In quiescence the accretion rate of our 
model (3.0 x lO^^Mg/yr) is similar to those of class I YSOs 
(10 '' - 10-^ Me/yr, White & HiUenbrand 200-fll or class II ob- 
jects (10 - IQ-^Me/yr, Gullbring et al. 11991 Hartmann et 
al. 1998). During the outburst it is 3.0 - 7.0 x IQ-^Ms/yr, 
higher than that of quiescent low-mass YSOs, comparable to 
that of FUors (10 - lQ'*Mo/yr, HK96). Intermediate values 
are expected for EXors (HK96). Juhasz et al. (I2012I I derived 
2.2 X lO^^Mo/yr for the prototype EXor, EX Lupi, at its out- 
burst peak in 2008. 

In the quiescent model we derived a similar accretion rate 
value as AP07 and Muzerolle et al. (I2005I I found, i.e., a few 
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times 10"^ Mo/yr. Aspin et al. (I2008I I claimed, based on the 
strength of the Bry emission line, that the accretion rate was 
somewhat higher IQ^^ M^/yr) after the outburst. However, 
the estimates of the accretion rate derived from the Bry mea- 
surements strongly depend on the applied value of Ay. Aspin et 
al. (12008 ) considered higher visual extinction than, e.g., AP07. 

As a final refinement to fit the NIR color variations, we also 
had to modify the parameters describing the accretion charac- 
teristics. At the outburst peak, the hot spot covers a large part 
(«i50%) of the stellar surface. In quiescence, the hot spot is much 
smaller, i.e., <1% surface fraction. Similar results were deter- 
mined for other eruptive YSOs (EX Lupi, Juhasz et al. 2012 
and V1118 Ori, Audard et al. i2010< . Calvet & Gullbring (.1998 j 
found fractions of 20% in strongly accreting systems, and typi- 
cally 1% for T Tau systems. These surface fractions correspond 
to 6500 and 15000K hot-spot temperatures at the outburst peak 
and quiescence, respectively. 

5.2. Dynamics of the inner disk and envelope 

During the outburst the temperatures were the highest at the 
peak, which defines the radius up to which dust had to evapo- 
rate. The variation of the inner radii of the disk and the envelope 
suggests that the physical mechanism enlarging the inner dust- 
free hole from 0.5 to 0.7 AU is the evaporation of dust by the 
increased heating from the central source. The intensified wind 
from the center of the system might have also played a role in 
these changes (Reipurth & Aspin |2004] |. 

5.3. Burst of luminosity or drop of extinction? 

Owing to the a structural change described previously, the ex- 
tinction of the system should also change. The increased ac- 
cretion luminosity and the varying extinction determine the to- 
tal brightening of the object. When the heat because of the in- 
creased accretion luminosity evaporates dust particles and cre- 
ates a larger inner hole, the line-of-sight extinction decreases, 
which leads to an additional apparent brightening of the object. 
Since the density in the envelope decreases with radius, most of 
the extinction comes from the innermost parts. Therefore, dust 
evaporation at the inner edge of the envelope can significantly 
change the total extinction toward the central star. Similarly, the 
gradual fading of the central source may allow re-condensation 
of dust grains in regions where evaporation took place previ- 
ously, resulting in an extra dimming. Kun et al. (1201 II) claimed 
that such re-condensation was observed at another YSO, PV Cep 

For the understanding of the outburst, it is essential to 
separate to which extent those two agents contributed to the 
brightening of VI 647 Ori throughout the eruption. Some au- 
thors (e.g., McGehee et al. 2004 Reipurth & Aspin I 2004l l con- 
cluded that only part of the brightness variation of VI 647 Ori 
was due to extinction variation, e.g., dust-clearing. From our 
models, we derived a change of the extinction AAy ~ 4.5 mag 
or AAij, X! 2.5 mag. Compared to the largest measured varia- 
tion in Ic (6.0 mag, AP07 and Aspin & Reipurth 2009 ), it leaves 
more than 3 mag intrinsic brightening in Ic. There are a num- 
ber of different attempts in the literature to derive the value of 
extinction and its variation during the outburst. Our model re- 
sult for the quiescent phase {Ay x!23.4mag) agrees quite well 
with the value derived by Aspin et al. (2008 Ay x 19 mag). 
The object moved along the reddening path in the NIR color- 
color diagram during the outburst. The excursion corresponds to 



AAy ~ 5 mag (e.g., Reipurth & Aspin l2004l l. which is close to 
our results. However, with different methods, Aspin et al. (120081 ) 
derived AAy ~ 10 mag. Since the optical extinction values are 
poorly defined over the outburst, we do not aim at fitting them 
with our model, only attempt a qualitative check of this param- 
eter. Our modeling indicated that the extinction is higher in qui- 
escence than at peak brightness. 



Table 3. Varied model parameters for different epochs. 



Parameters 


2004 


2005 


2003/06 




Mar 


Oct 


Mar 


Sept 


quiescent 


M (Moyr-' x 10^^) 


7.0 


5.5 


3.5 


1.6 


0.3 


/?in,disk (AU) 


0.7 


0.7 


0.7 


0.7 


0.5 


Rin.cn. (AU) 


0.7 


0.7 


0.7 


3.0 


0.5 


Ay (mag) 


18.9 


18.9 


18.9 


11.5 


23.4 



5.4. Scenarios for variations seen by MIDI 

A blown-up spherical cavity The expected structural changes 
can either occur relatively slowly during the March-September 
interval or in a few weeks' time (cf. the light curves in Fig. [TJ. 
Even in the latter case such a change of the envelope seems 
to be plausible (corresponding to 100 - 200km/s velocities). 
During most of this period the object was slowly fading, there- 
fore we would expect the decrease of the inner (dust evaporation) 
radii both of the disk and the envelope rather than any increase. 
Combet & Ferreira (12008) showed that at the accretion rate of 
10"'' Mo/yr, a strong disk wind, which can be the source of an 
outflow, is launched from the inner disk up to 2-3 AU. Thus the 
disk wind might be responsible for the clearing of the inner en- 
velope. However, during the 2005 March-September period the 
wind is expected to become weaker since its strength should be 
proportional to the accretion activity. Furthermore, the wind that 
was present from the beginning of the outburst and had no such 
effect before might only have cleared out the cone in the reflec- 
tion nebula when the outburst started (Reipurth & Aspin 20041 ). 
Therefore the stellar or disk wind is unlikely to have produced 
this clearing. 

Another possibility is that some temporary outflow changed 
the inner structure of the system and moved the inner edge of the 
envelope outward. Spectroscopic observations from early 2006 
showed evidence of a short outflow activity (Brittain et al. 2007X 
Furthermore, detection of forbidden optical lines indicated the 
presence of shocked gas in early 2006 (Fedele et al. i2007) . These 
lines can be tracers of a Herbig-Haro object (HHOfl If such 
an event were responsible for the change of the inner radius of 
the envelope (Rin.em), based on the beforementioned findings, it 
should be connected to the final decrease of the accretion rate. 
However, Brittain et al. (120071 ) proposed that because the inner 
radius of the accretion disk decreases together with the decreas- 
ing accretion rate, the re-structuring of the inner disk that pro- 
duced the (CO) outflow might have occured later. 

During an outflow activity or the appearance of an HHO, 
the disk structure can be very different than that of a standard 
accretion disk. There is a significant difference in all parameters 

Although HHOs are thought to be related to mass accretion, they 
are not typical in the environment of eruptive YSOs (HK96). Note that 
Eislijffel & Mundt ( 1997 ) identified V1647 Ori as the driving source of 
HH23 - presumably ejected some thousands of years ago - which is 
located 155"noi1:h of the young star and is close to axis of the nebula. 
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between those of a standard accretion disk and a jet-emitting 
disk (Combet & Ferreira 2008'(. Reproducing such events are 
beyond the limits of our static modeling and can be a source of 
discrepancy. 

By 2006 February, the object returned to its quiescent bright- 
ness at all wavelengths. This suggests that the structural changes 
in the circumstellar environment of VI 647 Ori should have been 
reversible on the timescale of some months. It means that the ma- 
terial should fill the evacuated spherical cavity on this timescale 
via accretion, which would require an inward radial velocity of 
10-20 km/s. In our modeling it means that the inner radii of both 
the disk and the envelope should move back close to the qui- 
escent values. Matter of the disk and the envelope might have 
moved inward due to the accretion and thus could fill the cleared 
regions. However, the dynamical timescale at ~3 AU is about 
seven years, much longer than half a year. Therefore dust re- 
condensation might also have played a role here. 

The disappearance of a warm halo In the next scenario we 
again consider a large spherical inner cavity in the envelope as 
in Sect.|43](/?in,env = 3.0 AU while /?in,disk = 0.7 AU), but we as- 
sume that the cavity has been produced at the early phases of the 
outburst for instance by stellar or disk winds (see e.g., Clarke 
et al. l2005l l. If dust was lifted from the disk above it by wind 
(Vinkovic & Jurkic|2007 Sitko et al. 2008 ), the dust could fill the 
inner cavity of the envelope. Note that such a halo cannot be dis- 
tinguished from the inner envelope in our models, but represent 
a different mechanism. The warm dust halo present in the system 
during the plateau phase makes VI 647 Ori look more compact at 
the first MIDI epoch. The apparent change of the visibility data 
could then be caused by the waning wind. As the accretion pro- 
cess and thus the wind becomes weaker, less and less dust is fed 
into the cavity, while the dust moved there previously is blown 
farther outward. In this case most of the dust - which made the 
object compact in 2005 March - is cooled below the tempera- 
tures corresponding to the mid-infrared wavelengths. By 2005 
September the cavity is practically cleared. Without carrying out 
new modeling, we can consider the results of the radiative trans- 
fer modeling described above. This scenario could explain the 
transition between 2005 March and September seen in the data 
and also the short timescale variability of the mid-infrared emis- 
sion. 



An out-of-equilibrium system Finally, we further speculate to 
find other alternatives for our series of static models. Between 
2005 March and September the heated disk area should shrink 
gradually because of the fading of the central illuminating 
source. In contrast, the visibilities of VI 647 Ori apparently did 
not follow such a scenario. The above hypothesis would pre- 
dict a smaller - thus less resolved - source at our second MIDI 
epoch. Our measurements, however, revealed the opposite be- 
havior: the visibilities decreased between the first and the second 
epochs. Since the shapes of the two visibility curves are simi- 
lar (Fig. 1.), and their ratio is almost independent of the wave- 
length (^1.4-1.5), a simple qualitative picture might explain the 
data. One could assume that the system consists of two compo- 
nents: one compact central and one extended component (e.g., 
approximated by a Gaussian brightness distribution as in Sect. 
3.1). If the resolved component did not change between the two 
MIDI epochs, but the emission of the central unresolved source 
had significantly dropped in the same period, the emission of 
the system would have become less peaked, i.e., relatively more 



resolved. This might be a simplistic picture of what caused the 
decrease of the measured visibilities. 

The constancy of the extended emission component could 
easily be explained if the ~ 10 yum flux arose from optically thick 
regions, whose temperature is not adjusted rapidly to the chang- 
ing central illumination field. Although we see variations of the 
mid-IR brightness of the system on weekly timescales, that of 
an optically thick component could not change significantly be- 
cause the corresponding timescales (Chiang & Goldreich 1 19971 1 
exceed even the difference between the epochs of the MIDI ob- 
servations. However, according to Muzerolle et al. (12005) . the 
contribution of an optically thick accretion disk to the mid- 
infrared emission of VI 647 Ori in outburst is small, and the 
~10//m flux mainly arises from an envelope. In this envelope 
the temperature of the dust grains must be quickly adjusted to 
the external radiation field. This suggestion is supported by the 
fact that the mid-infrared flux of VI 647 Ori had increased re- 
markably in less than a few months, between the beginning of 
the eruption and the Spitzer measurements in March 2004. Thus 
the invariability of the extended emitting component in the sim- 
ple picture above is not straightforward to explain. 

A qualitative argument leads to a possible solution in which 
the decrease seen in the MIDI visibilities is caused by the fading 
of the central source, which is not immediately followed by the 
fading of the outer regions. However, this effect should be con- 
firmed by dynamical modeling and is beyond the scope of the 
present paper. 

6. Conclusions 

We performed interferometric and photometric observations of 
VI 647 Ori during its 2003-2006 outburst to investigate the tem- 
poral evolution of its circumstellar structure and physical pro- 
cesses related to the eruption. In addition to the general fading 
of the object - shown by the multi-wavelength photometric data 
and archival mid-infrared spectroscopy, - short timescale varia- 
tions were also observed. Optical-infrared SEDs at five epochs 
were compiled. Our radiative transfer modeling, with a smoothly 
decreasing accretion rate as a major varying parameter, provided 
good fits of the SEDs at different stages of the outburst. It is im- 
portant to note that the inner radii of the dust disk and envelope 
also had to increase during the transition from quiescence to the 
outburst peak. The latter finding is clear evidence of dynamical 
variations in the inner circumstellar environment of V1647 Ori. 
This dust clearing is likely caused by the evaporation of the dust 
grains due to the outburst heat. 

High angular resolution spatial information were also con- 
sidered in the model-fitting procedure. VLTI/MIDI data obtained 
at two epochs, during both the slow and rapid fading stages, 
show a considerable change of the circumstellar structure. In 
constrast to our expectations, based on our model sequence, the 
object looked more resolved at the second epoch. One possible 
explanation can be a rapid removal of dust from the inner 3 AU 
of the envelope, possibly caused by wind or outflow processes. 
In one case this spherical cavity is produced at the end of the out- 
burst phase (the blown-up cavity scenario), or alternatively at the 
beginning (our warm dust halo scenario). Finally, we may also 
see VI 647 Ori in a non-equilibrium situation in 2005 September, 
when the sudden fading of the central source was not yet fol- 
lowed by the fading of the optically thick circumstellar material. 

In general, our modeling showed that the circumstellar envi- 
ronment of VI 647 Ori can be described by a disk and envelope 
system with parameters that are typical for embedded low-mass 
YSOs. This finding supports the hypothesis that eruptive YSOs 
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are not peculiar objects, but represent an important phase in the 
evolution of all low-mass YSOs. 
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Appendix A: Data reduction with different versions 
of IVIIA+EWS 

Since we aim at investigating temporal changes of the source 
structure, and the first MIDI datasets were reduced with the 1 .3 
version of MIA-nEWS, we repeated the data reduction with the 
1.6 version. It turned out that longward of 11/im the calibrated 
visibilities differ significantly with respect to the previous reduc- 
tion, i.e., are higher than were shown in Paper I. At shorter wave- 
lengths the results are consistent. Nonetheless, the results with 
EWS and MIA agreed well in Paper I. The EWS and MIA re- 
sults are also consistent for the second epoch data (Sect. lZTI ). As 
another test, we also reduced the second MIDI dataset with the 
older versions of the software. The results show a similar bias. In 
the case of MIA, the difference comes mainly from two changes 
mage since the first release. First, newer MIA versions use the 
photometry routine of EWS instead of its own (used in earlier 
versions). Second, these new versions of MIA fit fewer parame- 
ters of the mask than the older MIA version did. In the case of 
EWS, the difference is that its newer versions handle the data 
where the OPD difference is close to zero, differently from older 
EWS versions, i.e., these data are not considered anymore. Since 
all these changes of the softwares were made to improve the data 
reduction, we present calibrated visibility curves for both epochs 
produced wih MIA (1.6) and considered these for the radiative 
transfer modeUng throughout the paper. 
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